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Abstract 

0{a) QED radiative corrections to neutral current deep inelastic production of heavy 
quarks are calculated in the leading log approximation and compared with the corre- 
sponding corrections assuming a massless charm parton. Besides the inclusive case, cor- 
rections to the semi- inclusive d^a / dxdydz and the effect of z-cuts are studied. In the 
latter case, the massless corrections differ from the correct massive radiative corrections 
to deep inelastic heavy quark production by about 40%-10% for 0.2 < -2 < 0.5. 



1 Introduction 



Deep inelastic production of heavy quarks at HERA is important for several reasons: 

• At small X (x ~ 10^'^) the charm contribution F2 to amounts to 20%-30% 
making a reliable theoretical treatment of F2 necessary for a precision measurement 
of . 

• The bulk of heavy quarks is produced via the photon gluon fusion mechanism [|l], |^ 
providing the possibility to constrain the gluon distribution giy^jj?) in the proton 

ail- 

• By measuring differential distributions, the charm production mechanism can be 
tested, and it can be studied whether and when the charm quark behaves like a 
massless parton |]. 

In order to estimate radiative corrections to heavy quark production, we employ the 
well-known leading log approximation (LLA) 0, ||, . See also ref . [Jl0[ where radiative 
corrections to heavy flavour production have been studied in the Weiszacker- Williams 
approximation (WWA) using hadronic variables. The 0{q) (QED-)corrections to the 
process e + g — ^ e + c+c are shown in fig. |1] (not shown are crossed and virtual diagrams). 
The main part of the corrections comes from bremsstrahlung from the electron line (plus 
corresponding virtual corrections), see a) and b). These two contributions are enhanced 
by a large collinear logarithm Lg = \n{Q'^/m?^). In the same way, calculating c), d), and 
e), one gets logarithms Lc = \n{Q'^/'m'^) which are much smaller than Le because of the 
large charm mass m^, compared to the electron mass mg. For example, for = 10 GeV^, 
one finds Lg ^ 17.5 >> Lc ~ 1.5 {rric = 1.5 GeV). Diagrams d) and e) do not contribute 
for another reason: the photons are not resolved from the final state jets in general. (Only 
a light intermediary quark in diagram c) could give rise for a collinear singularity, which 
would have to be handled the same way as the corresponding collinear gluon emission. 
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leading to a negligible (at least in LLA) modification of the parton distributions.) Thus, 
we only take the leptonic corrections a) and b) (in 0{a)-LLA) into consideration. 

The paper is organized as follows: In Sect. 2 the formulae needed to calculate radiative 
corrections to neutral current heavy quark production in 0{a)-LLA are collected. In 
Sect. 3 numerical results are presented. Finally, the main results are briefly summarized 
in Sect. 4. 



2 Formalism 



We calculate the 0{a) QED-corrections to deep inelastic production of heavy quarks in 
leading log approximation (LLA) |^, ^. Following the formalism of ref. the cross 
section for deep inelastic production of charm quarks via the photon gluon fusion (PGF) 
mechanism (e + g^e + c + c) reads: 

da{ep ^ eccX) = / dzi / dz2 / dz3De/eizi,Q^)g{z2, fJ.'^)De/eiz3,Q^)dcr{e+g ^ e+c+c) 
Jo Jo Jo 



The functions -Dg/e, -Dg/e contain the factorized mass singularities from the electron line 
and are, up to 0{a), given by 



De/e{x, Q'') = De/e{x, Q') = 5(1 - x) + ( ^ ) PU^) 



with 



T -1 

Le = In — 



and the splitting function 

Peeix) 



1 



1 — X 

The (+)-distribution is defined by 

1 + z^ 

\ dz Pee{z) f{z) = / dz --— 
Jx J X ^ ^ 



(r3^ + 2^(^-^)- 



(2) 



(3) 



(4) 



(/(^)-/(l)) 







dz 

1 — z 



/(I) • (5) 
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De/e{zi^ Q"^) can be interpreted as the probability of finding an electron with momentum 
fraction zi inside the incoming electron, De/e{z3,Q^) describes the fragmentation of the 
scattered electron into the observed outgoing electron with momentum fraction z^, and 
g{z2,fi'^) is the gluon density inside the proton. Finally, da ist the cross section of the 
underlying hard scattering subprocess. 

The variables x, y and are reconstructed by the 4-momentum p'^ of the observed 
electron and the 4- momenta Pe, p of the initial state electron and proton: 

qi=Pe- Pe. Qi = -Qi , yi = , xi = = ^ , 6 

P-Pe ^p-qi Syi 

where S = {p + PeY- In the same way, one defines for the hard scattering process 
e{pe = zipe) + g{p = Z2p) — > e(]3g = p'^/ z^) + c{pc) + X: 

Q = Pe- p'e, = -q^ = —Qh y = ?^ = ^ — L±^^ (7) 

Z3 PPe Z1Z3 

X = 7-—: = = —f n V' S = iPe+P) = Z1Z2S . (8) 

2p-g Sy Z2[ziZ3 - I + yi) 

Because we only use leptonic variables x/, y/, etc. , the index / will be suppressed from 
now on. 

In the following, we calculate the radiative corrections to the photon gluon fusion 
(PGF) subprocess 

e(Pe) + g{p) — > e{p'^) + c(pc) + c{Pc)- (9) 

and compare the results with the radiative corrections to the charm excitation (CE) 
subprocess 

eiPe) + c(p) eip'^) + c{p,) (10) 

where, in contrast to the PGF, the charm quark is treated as an intrinsic (massless) parton 
of the proton. The latter process has been used by the HI Collab. ||rT| in a recent analysis 
of leptoproduction of D-mesons employing some charm density c(x, /i^). 
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2.1 Inclusive case 



2.1.1 Photon gluon fusion (PGF) 



The hard cross section for the process (j^ can be written as: 

(fa _ ina'^S 
dxdy Q4 



with the partonic structure functions WA [I3i 0, [15 



2£ 



(/2 - h) 



i-eicZiA.Q')OiW'-Ami), k = 2,L 



(12) 
(13) 



i_„(o) 
47r 



i_„(o) 
47r^'^ 



^ - a;2(l - x) + 2^x^(1 - 3x) - l^x^ 



In 



1 + /3 
1-/3 



4x (1 - x) 



X m 



2'^x''(l-x) 



1-p 



2 Q2 

+ 2/?x2(l -x), 



(14) 
(15) 



where Cc, mc are the charm quark charge and mass, respectively, ^ = ^(Q^ 



m 



2 ' 



W\x,Q' 



Q\l-x) 



X 



and = /32(x, Q^) = 1 - Im^/iy^ 



By insertion of eq. (0), (0) into eq. (|T|), one obtains the cross section da = da^ + 
da^ + da^ . The Born cross section da^ is given by 
c/VO Ana's 



dxdy Q'^ 
with the leading order structure functions 

dz. 



[{l-y)F,{x,Q')+xy'F,{x,Q')] 



J ax ^2 

F2(x,(5^/i^m^) = /" —Z2g{z2,n'^)f2ix/z2,Q',ml,i/), 

J ax ^2 



(16) 

(17) 
(18) 



where a = 1 + Arn^/Q'^. da^'-^ contain the contribution due to radiation of a colhnear 
photon from the incoming (initial state radiation, ISR) and outgoing electron (final state 
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radiation, FSR), respectively: 



dxdy 



^ aL, 



dzi 



l + zl 



2% 



(ao(zi,l)-ao(l,l)) 



(19) 



dxdy 



+ 



27r 
aLe 



2tt 



dz:K 



1 + 4 
1-Zs 

o^o(l,l), 



(ao(l,Z3)-^o(l,l)) 



(20) 



with 



0-0(2:1,^3) 
H(z) 



dz2 g {z2, /i^) 
dx 

1 - a; 



d{x,y) d'^a 



d{x, y) dxdy 
21n(l 



z) +z+-z^ 



(21) 
(22) 



and the Jacobian 



y 



d{x,y) ^ 

d{x,y) Z2Z3{ziZ3 - I + y)' 



The electromagnetic coupling is taken to be a = a^Q"^ 
tion bounds 



(23) 

1/137.036[|. The integra- 



y 



+ 



Arrii 



z^^"'{zi, Z3) 



1 — xy S{1 — xy) ' 

1-y + xy 

1-Aml/S' 

Amlzs^ 
1 + I zixy- 



(24) 



1 



zij ziZ3 + y - 1 

follow from the conditions W"^ > 4m^ , < z'^'^^{zi, z^) < 1. If no photon is radiated 
{zi = 1, Z3 = 1), one finds the well-known expression z™"(l, 1) = (1 + 4m'^/Q'^)x = ax. 
In the kinematical region of small x, the bounds read approximately: 
1 — y + 0{x) + 0{m1/S). It should be remarked that ctq is related to the Born cross 
section by d'^a^ /dxdy = cro(l, 1). 

^ Using a running coupling a{Q^) according to eq. (14) in taking effective quark masses |^ 
m« = md = 0.041, rus = .15, rric = 1.5, and mb = 4.5 GeV, instead of a constant a, leads to negligible 
differences. 
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2.1.2 Charm excitation (CE) 



In the case of electron quark scattering, the radiative corrections in LLA are known up to 
|, 0, |, g lig. For comparison with eq. (||), (|0D, we use equations (22)-(29) in 0, 
with the charm quark as the only massless parton (MP) in the initial state. Furthermore, 
contributions due to Z-exchange can be neglected in the realm of < 50 GeV^, i. e. , 
we make the replacements — e^, — » (in eq. (22)-(29) in [0). For the partonic 
structure functions, the simple relations 

f2{x,Q') = el6{l-x) , (25) 

W,Q') = ^f2{x,Q') (26) 
2x 

hold, and one easily finds [0 

ao(zi, ^3) = ^^[1 + (1 - y f]el Uz-2, Q') + -c{z2, Q')) (27) 
ziz^y-^b ^ ' 

with 

-Z2 = ^i^^, ^— = 1^, = 1-^(1 -a;) . (28) 

^1^3 + 1/ - 1 \-xy 

(Of course, eq. (p7|) and ( pH]) have to be inserted into eq. ([T9|) and (pO]).) 

Finally, we briefly discuss the Compton contribution da'" to the radiative corrections 
which can be obtained from eq. (21) (Qo = 200 MeV) and (29) in 0, only allowing 
for the charm (or anti-charm) quark in the initial state. Fig. displays the Compton 

contribution 5*" = - — — ;~ fo^^ x = 10^^ (dashed line) and x = 10^^ (full line) using 

dxdy dxdy 

the CTEQ4L parton distributions il^- For y < 0.7, 5*" is small and can be safely neglected 



for experimentally relevant values of ?/ 
2.2 z-difFerential case 

In order to calculate radiative corrections to the ^-differential cross section d^a^ / dxdydz, 
where z = p ■ po/p ■ q, the fragmentation of the charm quark into the observed Z)-meson 
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must be taken into consideration. Thus, it is necessary to extend eq. (|1|): 

da{ep^eDX)= / dzi / dz2 / dz3De/e{zi,Q'^)g{z2, iJ.'^)De/e{z3,Q'^) / dz^Ddz^jda. 
Jo Jo Jo Jo 

(29) 

The hadronization of the outgoing charm quark with momentum pc into the observed 
D-meson with momentum p£, = z^p^. is modeled by the fragmentation function Dc{zi). 
Writing the partonic cross section differential in Zc = p ■ Pc/p ■ Q, the relations in eq. (|^) 
have to be accomplished with 



zyz3 



_ z 



with r defined by 



Z4{ZiZ3 + y-l} Zi 



y Z3 



Z1Z3 + y-l 

which can be derived from the definitions oi Zc, z and g, using po = Z4Pc. 



(30) 



(31) 



In analogy to the above discussion, we compare the "massive" (PGF) with the "mass- 
less" (CE) corrections, using a massless charm parton (MP). 

2.2.1 Photon gluon fusion 



The partonic cross section reads 
d^a ATia^S 



1 - y)f2 ( X, Q^; z^ ) + xy^fi ( x, Q'^; z, 



dxdydzc Q4 
where the structure functions are given by ^ 

'1 



/2(x,Q2;C) 
/l(x,Q^C) 



«s 2 



^BG (x, Q\ C, m,) + -BL (x, Q^ C, m,) 



TT Z 



with 



BL (a;,Q^C,"^c) = 4x 



1 — X — X- 



rrir 



Q' C(i-C) 



^The corresponding expressions for general boson gluon fusion can be found in ||H 



(32) 



(33) 
(34) 



BG (x,Q^C,"^c) 



.rrir. 



l-2x + 2x^ + 4— - x) 

Q 



C(i-C) 



■\-2x^-^ (1-2 



Q2 j (1_C)2^2^ 



ra, |T8[ and Q ; ic) = — /2(x, Q ; 4) - Q \ ^c) 



2x 



The cross sections for ISR and FSR have the same structure as in eq. (|T9D and (pO|) 
with a z dependent function (Tq: 



dxdydz 
dxdydz 



dzi 



l + zl 
1-^1 



27r , 

+ ^iJ(^r>o(i,i;^) 

ZTT 



((To(^i,l;2) -cro(l,l;^)) 



(35) 



dz'x 



27r , 



1+i 

1-^3 



{aQ{\,zj,;z) -cro(l,l;z)) 



(36) 



with 



0-0(^1,^3; z) 



dzS^M^J^D^^zi) . (37) 



9(x, z) dxdydzc 



The Jacobian in eq. ( P7D can easily be calculated as 

d{x,y,Zc) d{x,y)dzc 



r 



d{x,y,z) d{x,y) dz ^2^4(^12:3-1 + 2/)^ 



(3^ 



The integration bounds can be derived from the conditions W"^ = ^ > 4m^, lead- 



ing to eq. (0), and z™" < < ^c"'"'' with zi 
1 — 4ml /W"^, implying 



l±/3(x,g2 



, where /3^(x, 



^r"" = min 



mm 



zr 



ZT 

max 

^ ' Zmax 



(39) 
(40) 



with r from eq. (pT]). For Dc{z), we use a Peterson et al. fragmentation function ||20| 

D,(^) = iv|^[l_^-i_£^/(l-^)]2|"' , (41) 



normalized to dzDc{z) = 1, i. e. 



(£,2_6£, + 4) 



(4 - ec)V4i7^i7 



arctan ■ 



VAe. 



c '~c 



1 1, 1 

+ arctan = > H — m H 



(42) 



2.2.2 Flavor excitation 

In the subprocess e + c — > e + c we have Zc = p ■ Pc/p ■ <? = 1, following from Pc = ^ + P- 
Thus, one obtains the Zc-differential cross section by inserting the structure functions 



/2(x,Q',^c) = elS{l~x)S{l-z,) 
into eq. Because of 5(1 — Zc) = Z46{z4 — rz) and dzc/dz = r / z^ one finds 



0-0(21, 23; z) = aoizi, Z3) rD^{rz) 



(43) 
(44) 



(45) 



with cro(zi, 2:3) from eq. (|27| ) and r from eq. (pl|). 



Finally the z-differential corrections are given by eq. (p5D and (|36| ) with (Jq from eq. (|45| ) 
and the integration bounds 



max 



1 — 



,1-1/(1-2:) 



max 



1 - 1/(1 - a;), 



l-y 



The second boundary in max[ 



1 - yz_ 

can be deduced from < rz < 1. 



(46) 
(47) 



Eq. (|35|), (|36|) can easily be tested: 



dz- 



(48) 



/o dxdydz dxdy 

with d'^a^'f /dxdy from (|T^, (PU]). The Compton contribution will be neglected for the 
same reasons as in the inclusive case. 
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3 Numerical results 



As usual, the radiative corrections will be shown in form of a correction factor 6 defined 
by da = da°{l + 6), i. e., 6 = 6' + 5^ + . . . with 

' dxdy dxdy 

or, in the 2;-differential case, 

5''^{x,y,z) = , , , /- 



dxdydz dxdydz 

In all figures we use HERA centre-of-mass energies S" = 4 ■ 27.5 ■ 820 GeV^. 

Fig. 1^ shows the radiative corrections to heavy quark production (PGF) for experi- 
mentally relevant values^ of Bjorken-x as a function of y according to eq. (|19D and (|20|) 
using the GRV94(LO) parton distributions |PT[. The factorization scale has been chosen 



to be fi^ = + Ami lH with = 1.5 GeV. For x = 10~^, the typical shape of radiative 
corrections in leptonic variables can be seen, with large corrections for ?/ ^ 1, whereas 
for smaller x < 10~^ the curves become more and more flat for y ^ 1. 

The theoretical uncertainties due to different choices of parton distributions, factor- 
ization scales, and charm masses turn out to be small, as can be seen in fig. ^ and fig. |^. 
One finds 5(CTEQ) - 5(GRV) < 0.02, S{mc = 1.3) - 5{mc = 1.7) < 0.03 for relevant 
y < 0.7, and = Ami) ~ "^(/^^ = Q"^ + Ami) < 0-03, where the scale /x^ = Ami has 
been favored in This could be expected because variations of /i, mc, or the parton 
distributions lead to rather similar changes in da^ and da"^'^ , so that the quotient 5 does 
not change too much. 

In the recent analysis of deep inelastic charm production by the HI Collab. ||l[, the 



radiative corrections have been calculated in 0{a)-LLK with the HECTOR package 



using the charm excitation subprocess (F£ = 0) and the GRV92 parton distributions . 



Furthermore, only initial state radiation has been taken into consideration, assuming the 
•^For a clean extraction of the gluon density it is necessary to extend the x-range to x < 5 • 10^'' |^. 
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collinear final state photon not to be separated from tlie outgoing electron |TT|, leading 
to small corrections. Because PGF has been measured to be the dominant (> 95%) 
charm production mechanism in the small x (and Q"^) range it is necessary to check 
if the "massless" corrections (5*^'^) agree with the "massive" ones (5^*^^). In fig. ^ we 
compare the massive (/i^ = + 4m^; =1.5 GeV) with the massless corrections due 
to initial state radiation. The experimentally relevant values of ^ are indicated by 
dotted vertical lines. The conversion of "masslessly corrected" (5*^'^) data to "massively 
corrected" (5^'^^) data can be performed by applying the factor R = (1 + 5*^^)/ (1 + 5^*^^) 
because of 

da{PGF) = da{exp) ^ 



1 + 

1 + 5^P 



= da{MF)^^^^ = da{MF)R. 

For the (x, Q^) data points one finds |i? - 1| = - 6^^^\/{l + 6^^^) < 3%,i. e. , the 
differences between massless and massive radiative corrections lead to a small increase 
^^MP > (J-P"^^) of the data. However, this difference is small enough to use the simpler 
charm excitation subprocess for calculating the radiative corrections in the inclusive case. 

Of course, heavy quark production processes are exclusive in the heavy quark momen- 
tum and on this more differential level the photon gluon fusion and the charm excitation 
processes are not compatible which can be seen, e.g., from the different shapes of the 
xd = IPdI/IPpI = 2\'p^\/W distributions shown in the experimental analyses ^ || (fig. 6, 
fig. 1 resp.). We prefer to employ the lorentz-invariant variable z = p ■ po/p ■ Q which 
approximately transforms into Xd in the 7*-p centre-of-mass system, following the argu- 
mentation in |[T|. This can be easily seen by calculating z in the 7*-p-CMS. One finds 
z = XDsm^9*/2 with]] 9* = ^{p*ph) ~ tt. 



In fig. we show the z-differential radiative corrections according to eq. (p5D, (|36D, 
]) (PGF, full line) and (|45|) (MP, dashed line), where it has been integrated over the 



+Pp = Oi Pg.c — ^2Pp ^ ~^Pg,c = ~(1 ~ ^2)Pp- This means for the CE-subprocess: p* + p* — 
p'l = Z4p5 ^ 0* ^ -k and for the PGF: p* + p* = p* + p^ ^ 9* ^ n (the coUinear configuration is 
dominant). Note that an angle 6** = 160° stih yields sm^e*/2 = 0.97. 
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kinematical range 10 GeV^ < < 100 GeV^ 0.01 <y < 0.7. The left hand side shows 
the correction factor 5* = da^/ da^ as well as the sum of initial and final state radiation 5* + 
5-^. In all cases, we employed the GRV92(LO) (massless) and GRV94(LO) (massive) 
pT| parton distributions, in the massive case with rric = 1.5 GeV and ^'^ = + 4m^. 



For Dc a Peterson et al. fragmentation function ||20[ with e = 0.15 has been taken. To 
study the dependence of 6 on e, we compare the radiative corrections for two different 
choices of e in fig. |. In the massive case the e-dependence is obviously small, whereas 
for massless corrections one finds big differences in the steep region z < 0.4. Fig. ^ 
reveals a rather big difference between massive and massless corrections for z < 0.5. 
For z —>■ 1 only soft photon radiation is allowed, so that the corrections factorize and 
become independent of the underlying subprocess. On this more differential level it is 
necessary to calculate the radiative corrections using the photon gluon fusion subprocess 
because of the big deviations of the massless from the massive corrections for z < 0.5. 
On the right side of fig. |^ the 2;-dependence of the cross sections da^'^'-f is displayed. 
Integration over < z < 1 leads to the observed small differences between massless 
and massive corrections because positive and negative contributions compensate each 
other. To get an impression of the effect of 2;-cuts, we show in fig. |^ the correction 



factor y] 2;™", 2"""^) = J^„i„ dz^^/ J^„i„ dz^^^ for three 2;-integration ranges, 

from top to bottom 2™" = < z < z'""^ = 1, 0.2 < z < 0.8, and 0.3 < z < 0.9 
for X = 10^^. For completeness, we used the GRV94/92(LO) parton distributions El 



| 23| for the massive {rric = 1.5 GeV, /i^ = + 4m^) and massless corrections and a 
Peterson et al. fragmentation function with e = 0.15. As can be seen from fig. |^, 
a cut z > 0.2 leads to the exclusion of positive contributions of rfo"*'-^, diminishing the 
radiative corrections. 



4 Summary 

The (9(a)-QED corrections to inclusive and 2;-differential deep inelastic electroproduction 
of heavy quarks have been calculated in the leading log approximation, using electron 
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variables. The results have been compared to the radiative corrections in the MP scheme, 
where the charm quark is assumed to be a massless parton in the proton. In the inclusive 
case, the differences l:)etween these two approaches turned out to be negligible. However, 
the measurement of heavy quark production is of course differential in the momentum 
of the (observed) heavy quark, recommending to perform the radiative corrections on 
the same differential level. Thus, we have considered the semi-inclusive ^^-differential 
case, in which the massive corrections have to be applied, i.e., using the photon gluon 
fusion subprocess, because for 0.2 < 2; < 0.5 the massless corrections differ from the 
massive ones by about ~ 40%-10%. Furthermore, we studied the effect of cuts on the 
2;-integration range. A cut z > 0.2, e.g., excludes positive contributions of da^'^ , so that 
the (2;"**" < z < 2;"*"^)-integrated corrections are smaller as in the fully inclusive, i.e., 
(0 < 2; < l)-integrated case. 
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Figure 2: Compton contribution 5'~^ according to eq. (29) in [Qq = 200 MeV, considering 
only the (anti-)charm quark in the initial state) with x = 10^^ (full line) and x = 10~^ (dashed 
line), employing the CTEQ4L parton densities [|17[. 
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Figure 3: Radiative corrections to heavy quark production (PGF) in (9(a)-LLA, using the 
GRV94(LO) parton distributions and the factorization scale ji^ = + 4m^ with rric = 
1.5 GeV. ("initial" and "final" refers to initial and final state radiation.) 
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Figure 5: Left side: 5{x,y), for two different factorization scales fi^ = + 4m^ (full 
line) and /i^ = 4m^ [|2^ (dashed line) {rric = 1.5 GeV). Right side: S(x,y), using the charm 
masses = 1.3, 1.5, 1.7 GeV (/i^ = 4m2).(PDF: GRV94(LO) 0].) 
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Figure 6: Comparison of massive (solid line, /i^ = + 4m^, = 1.5 GeV) and massless 
(dashed line) radiative corrections S (initial) for experimentally relevant x |]T| using the parton 
distributions GRV92(LO) (MP) O and GRV94(LO) (PGF) The dotted vertical lines 
indicate values of constant = 12, 25 or 45 GeV^. 
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Figure 7: z-differential radiative corrections to massive (rric = 1.5 GeV, /i^ = + 4m^) 
(full line) and massless (dashed) charm production, integrated over the kinematical range 
10 GeV^ < < 100 GeV^ 0.01 < y < 0.7, using the parton distributions GRV92(LO) 
(massless) |^ and GRV94(LO) (massive) and a Peterson et al. fragmentation function 
0] with e = 0.15. 
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Figure 8: As fig. ^ (Ihs) with e = 0.15 (full line) and e = 0.06 (dashed line). The contribution 
due to initial state radiation is shown (Ihs) as well as the sum of initial and final state radiation 
(rhs). 
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Figure 9: 5{x, y; 2™", ^max^ (defined in the text) in dependence of y for x = 10 ^ and three 
choices of 2™", z™"^: From top to bottom < z < 1, 0.2 < z < 0.8 and 0.3 < z < 0.9. The 
massive corrections have been calculated employing the GRV94(LO) parton distributions 



and the factorization scale = + Ami ^'^'^ = 1-5 GeV. For the massless corrections 
we have taken the GRV92(LO) parton distributions In all cases we employed a Peterson 
et al. fragmentation function with e = 0.15. 
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